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ABSTRACT 



Many (perhaps most) service systems, such as repair and job shops, 
computation centers, and transportation networks, experience demand that 
is non-stationary in time. This paper describes models for situations 
in which demands made are by a finite number of individuals, who, having 
been served, do not return until much later. Such a transitory demand 
or arrival process describes many phenomena, among them being commuter 
rush hours, and also perhaps the effect on a population of individuals 
their simultaneous exposure to a dosage of medicine, a disease, or even 
a pollutant. Our paper formulates several models for the service of such 
demands and describes the manner in which system state may be approximated 
by Gaussian processes, in particular the Orns tein-Uhlenbeck and Wiener 
diffusions . 
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I . Introduction , 

The theory of service systems and waiting lines has developed 
extensively over a number of years, stimulated by a desire to describe 
congestion and traffic in many applied situations. Almost without 
exception, however, the mathematical theory has dealt with stationary 
models, in which a steady stochastic flow of arrivals approaches a set 
of servicing facilities; for recent exceptions see the work of Newell 
[1(3, and that of Milch [9]. Nevertheless, many quite reasonable situa- 
tions suggest the study of models that are non-stationary . In this 
paper we describe and analyze a class of models of a distinctly non- 
stationary nature. We call them transitory service systems. 

Transitory phenomena occur when members of a finite population 
of customers makes application for service once, or at most a finite 
number of times. Some examples follow. 

(a) A questionnaire, letter, or newspaper or television advertise- 
ment offering a prize or job, is transmitted. A certain number of 
individuals choose to respond, and do so after various times. The 
responses must then be processed. Of interest is the number of indi- 
viduals in various stages: e.g. potential responders waiting to respond, 
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and responders undergoing processing. One may also wish to infer the 
number of potential responders in the population on the basis of the 
responses during a short initial time period. 

(b) A message is sent by radio to a group of stations or ships. 

The message requires a reply or verification, an act that takes place 
after a random time, and requires a random time to complete. The 
probability distribution of the number of messages being verified at a 
particular moment after transmission may be useful, as are other prop- 
erties of the reply process, 

(c) Certain copies of a new model of automobile or other consumer 
(or industrial or military) product contain a defect which is possibly 
dangerous and should be modified. When the defect becomes evident a 
recall or correction notice is sent out, following which the copies 
straggle in for modification. We are interested in the number of copies 
that have applied for, and are undergoing, modification, at any time. 

We are also interested in the length of time to complete all modifica- 
tions . 

(d) A ship puts out to sea with all systems in operative condition. 
Realistically, however, each system has a chance of failing during a 
voyage. Let us assume, as may hopefully often be reasonable, that at 
most one failure per system will occur during a voyage; it is this 
assumption that makes our process transitory. When the failure occurs 

a repair process is initiated. We are interested in total ship unavail- 
ability (the number of systems on repair) as a function of time, and 
other related measures of total system degradation. 
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(e) It has been suggested [2] that a road or freeway network may 
be considered to be analogous to an infinite server queueing system, 
with drivers f times on the freeway identified to be service times. 
Clearly, there will be certain intervals (rush-hour periods) during 
which a large segment of the population of automobiles will enter the 
freeway. They will not return until the next rush hour. Consequently, 
freeway occupancy appears to be a succession of transitory problems, 
one for each rush-hour period. 

Other examples will occur to the reader. We shall now embark 
on an analysis of the stochastic models suggested. Our particular 
emphasis will be upon approximations of a simple type. 
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II. " Sufficient Server" Models . 

Let there be N potential applicants for service in a population. 

The time of application of customer i, measured from some initial 

instant, is denoted by T.. The service time of customer i is S.. 

1 1 

We assume that at application each customer may immediately begin 
service, i.e. there are "sufficient servers." This is analogous to 
the infinite server model of conventional theory. 





Clearly applicant i 


is. 


at time t. 


a) 


not yet in service if 


T. 

l 


> t, 


b) 


in service if 


T. 

l 


< t, but T. + S. > t, 
9 i i 9 


c) 


finished with service 


if 


T . + S , < t. 
l i 



Consider the following models. 

Model 1 . {T_,} and {S_^} independently and identically distributed 

(i.i.d.) . 

Let T_^ have the distribution F, and the distribution 

G; let F(t) = 1 - F(t) , and G(t) = 1 - G(t) . Then if A(t) denotes 
the number of applicants who have not yet applied, Q(t) denotes those 
who have applied and are in service, and C(t) denotes the number of 
those who are finished at t we can write down directly, on the basis 
of independence, the joint generating function 

E[u A(t) w C(t) ] = (uF(t) + vG*F(t) + wG*F(t)} N . (2.1) 



where 



G*F(t) 



t 

* 

[1 - G(t-x) ]dF(x) 



( 2 . 2 ) 



0 



G*F(t) = 



G(t-x) dF (x) . 



(2.3) 



0 

If we put u = w = 1 we see that Q(t), the number in service has 
the generating function 

E[ v Q(t) ] = {[1 - G*F(t)] + vG*F(t)} N (2.4) 



of the binomial distribution. Therefore the expectation and variance 
of the number in service at time t are 



E[Q(t) ] = N 



[1 - G(t-x) ] dF (x) 



(2.5) 



Var [Q ( t) ] = N 



0 



[1 - G( t-x) dF(x) { 1 



[1 - G (t-x) ] dF (x) } . 

0 



As N becomes large one can approximate the distribution of Q(t) by 
the normal with mean and variance matching (2.5). 

The conventional infinite server problem, to which this model 
compares, yields no such simple time-dependent result. Recall that if 
arrivals are stationary Poisson then the stationary distribution is 
also Poisson; see Riordan [12]. Of course, if population size, N, 
is large then Q(t) may turn out to be approximately Poisson at a 
fixed time t. Consider the following extensions. 
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Model 2 . {T_^} and { S _^ } i.i.d.; N is Poisson. 

If population size N is a random variable with the Poisson 
distribution of mean a, then clearly Q(t) is also Poisson with 
mean aG*F(t) , as follows from the generating function: 



E[v Q(t) 



00 _ n 

] = l {1 - G*F(t) + vG*F(t)} n e 3 — 
n=0 n - 



= exp{aG*F(t) (1-v) } . 



( 2 . 6 ) 



Moreover, A(t) , Q(t) , and C(t) are jointly and independently 
Poisson in this model as is evident from the fact that the generating 
function factors: 



E[u A (t)vQ(t)wC(t) j _ eX p{ a p( t ) (i- u ) }exp{aG*F(t) (1-v) } 

exp{aG*F(t) (1-w) } . 



(2.7) 



It cannot, however, be concluded that {Q(t),t^0} is a Poisson process . 



Model 3 . {1\} are i.i.d. Gamma; { S ^ } are i.i.d.; the population 

size N becomes large. 

Suppose the arrival distribution, that of T\, is Gamma, with 

density 



f(t) = e 



-Xt (Xt )^ 1 

r(k) 



( 2 . 8 ) 



for t ^ 0 and k > 1. Then we can express the generating function 
of Q(t) , (2.4) , as 
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E[v Q(t) ] = [1 - (1-v) 



t 

- 

0 



,k-l 



n \ n -Ax N(Ax) , n N 

u ' G<t - x)lc irTio dxl 



(2.9) 



Now scale time by letting A N = 6, a constant, and allow N 
Then 

t 



E[v^^] -> exp{-(l-v)0 



k-1 

r - - , N 1 X dX > 

[1 G(t x)] r ^ k j } 



( 2 . 10 ) 



The continuity theorem then implies that under the stated 

conditions the distribution of Q(t) approaches the Poisson with 
-t k-1 



paramters 0 



[1 - G(t-x) ] 



dx 



0 



r(k) • 



Notice that if k = 1 then in 



the limit we have the familiar "infinite server" solution of the classi- 
cal theory: if arrivals are a stationary Poisson process — to which 

the process {A(t),t^0} converges as N « — then {Q(t),t^0} is a 

't 

[l-G(x) ]dx -> 0E[S] 



0 



non-s tat ionary Poisson process with E[Q(t)] = 0 

as t -* <». For other k values the process behaves in a comparable 

fashion. It may be contended that the type of large-N result just 

derived actually justifies the use of the classical queueing models 

that assume Poisson arrivals with stationary independent increments, 

at least early in the arrival process. Of course, an approximation 

based upon NX = 0 cannot be accurate late in the process, e.g. after 
Jin 2 



a time t^ 



-, when on the average one-half of the calling 



'1/2 A 

population has arrived. Therefore, in the following section we develop 
a procedure for approximating our processes A(t), Q(t), and C(t) by 

Ornstein-Uhlenbeck processes; these approximations can be expected to 
be reasonably accurate for moderate N, and t-values of practical 



interest . 
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III. Diffusion Approximation for Sufficient Servers . 

Since A(t) , the number of arrivals in time t, and Q(t), 
the number in service, are approximately normally distributed for large 
N, we are led to consider the use of a diffusion approximation for 
these processes; see Feller [4], and Cox and Miller [3]. We remind the 
reader that a diffusion is a Gaussian process with continuous sample 
paths. We next proceed to develop approximating diffusions for certain 
transitory systems, and then to make use of such approximations to 
obtain properties of the service systems. Our approach here is in the 
spirit of McNeil [8]; other authors have utilized similar procedures. 
For example, see Anderson and Darling [1], and Whittle [14]. 



III-l. Approximating the Arrival Process. 

Let A^(t) denote the number of arrivals by t out of a popula- 
tion of N. Each arrival time is independently distributed with d.f. 

F(t) , density function f(t), and hazard rate A(t) = f (t) [ 1— F ( t) ] 

V‘> 



Now clearly 



■+• F(t) in probability as N -*■ ® for every 



fixed t. In fact, the convergence is uniform with probability one by 
the Glivenko-Cantelli theorem; see Tucker [13], p. 127. Hence, we are 
led first to approximate A^(t) by a deterministic component , x(t), 
for large N: 

^(t) ~Nx(t) (3.1) 

Noting that if A^(t) arrivals have occurred by t, and thus that the 
probability of another arrival in (t,t+dt) is essentially 
[N - A^(t)]X(t)dt we can derive 



^ = X (t) [1 - x(t)], 



(3.2) 



the solution of which is x(t) = F(t). 
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To include a stochastic element, or noise , study 

A^t) - Nx(t) 



V c) ■ 






(3.3) 



Our first approach is to derive the transformed version of the forward 
differential equation satisfied by the density function of , pass 

to the limit as N -* <», and thus show that if converges to a 

limiting process the latter must be non-s tat ionary Ornstein-Uhlenbeck; 
see Cox and Miller [3]. 

If 

* 0^(0 _ 

(3.4) 



i 0 V t} r~ 

<j> N (0,t) = E[e ] i = /-l. 



is the characteristic function (ch.f.) of A^(t) we h ave 
i0Aft+dt) i6\,(t) 

E[e ] = E [ e (1 - [N-A^(t) ]X(t)dt}] 



+ E[e 



10 (Aj^O+l) 



{N-A N (t)}X(t)dt] + o(dt), (3.5) 



and upon collecting terms and letting dt -* 0 there emerges the 
differential equation 

9 V° i0 3 * n 

— — — = (e -l){N<)) N (0,t) + i -gQ-)X(t) 



(3.6) 



Now 



to derive an equation for the ch.f. of S^(t) , 



ieyt) 

+ N (t) = E[e ], 



(3.7) 



we note that 



<i) *„(8,t) - e iex(t) '^ *„ 



(ii) ^f(e//N,t) - e iex(t)^- 'V*' 11 + 16x .(t)/» e lex(t),¥ * N (8.t), 
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(iii) ^|(6/^,t) . e 19x < t>,/if ^|(e,t> + (3.8) 



Next change 0 to 0' = 0/t4T in (3.5), utilize (3.7), expand the 
exponentials - 1 = — - ~~ + 0( -jrr) ) and collect terms. The 

M 2N n 3/2 

essential outcome is 



9^. 



3 f + * N ie*¥{x'(t) - x(t) [l-x(t) ] } = - x(t)[i-x(t)]^ N 



1 

- ex(t) -rf + o(— ) . (3.9) 

90 v¥ 

In order for the equation to be satisfied by a limiting ch.f., tjj(0,t), 
as N -* 00 the term in brackets must vanish; solution of the resulting 
differential equation once again presents the deterministic equation 
of (3.2). There remains the equation 

H = - ex(t) If (3.30) 



which is the transformed version of the forward equation for a non- 
stationary Ornstein-Uhlenbeck (O.U.) process; see Cox and Miller [3], 
pp. 218-219. The infinitesimal mean and variance are respectively 

3(x,t) = -xA(t) and a(x,t) = f(t). Since the O.U. process is Gaussian 

0 2 

we may differentiate the Gaussian ch.f. exp[i0p(t) — — a 2 (t)] and 
equate coefficients in (3.9) to obtain the mean and variance. Not 
surprisingly. 



E[S(t)] = 0, and Var[S(t)] = F(t)[l-F(t)] (3.11) 



Hence we are led to approximate 



V fc > 



by 



Ajj(t) « NF(t) + /N S(t) 



(3.12) 



where S(t) is an O.U. process with parameters given by (3.11). 
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An alternative derivation of the O.U. parameters is easy, once 
one settles upon a diffusion approximation. The stochastic differential 
equation for A^ is given approximately by 

dAjj(t) = X(t)[N-A N (t)]dt + Z(t)/X(t)[N-A N (t)]dt, (3.13) 

Z(t) being a M purely random process," or the "derivative" of a Wiener 
process. The first right-hand term of (3.12) is the deterministic com- 
ponent or drift, and the second is the supplementary randomness or noise, 
the scale of which follows from considering possible events in a dt 
interval. Now apply normalization (3.3) to obtain 



dS N + v ^ x, ( t ) “ X(t){l-x(t)}] = - X(t)S N (t)dt 



+ Z(t) 



/ 



and once again we are compelled to set x T (t 
realize that the limiting noise satisfies 



A(t) 



i-x(t) - 



) = X (t) [l-x(t) ] 



dt, (3.14) 
and to 



dS (t) = - X(t)S (t) + Z(t)/f (t)dt, (3.15) 



the stochastic differential equation for the O.U. process with drift 
and variance parameters already derived. 



III-2. Passage Probabilities for the Arrival Process. 

In this section we show how the limiting O.U. process furnishes 
an approximation to the probability that A^(t) ^ B(t) for all t, B(t) 
being a suitably selected boundary. Our asymptotic result is closely 
related to the one-sided Kolmogorov-Smirnov test, and may be used for 
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the same purpose, i.e. as a test of the statistical hypothesis that 
the arrival distribution F(t) generates arrivals more rapidly than 
does F^Ct), i.e. that F(t) ^ F^(t) for a H t. The methodology 
employed may also be used to estimate the maximum queue size during a 
single rush period, but further approximations are required and this 
problem will be attacked later. 

It is shown in [3] that an O.U. process, S(t) may be expressed 
in terms of a Wiener process, X(t) : 

S(t) = a(t)X{x(t) } (3.16) 



where r(t) represents a time scale transformation, and a(t) is a 
real differentiable function. In order to represent our particular 
process we have 



a 1 (t) 
a(t) 



X(t), and a 2 (t)x'(t) = f(t). 



(3.17) 



Solving, we find that 

a (t) = 1 - F(t) , and x(t) = ^ -^F( ~ t ' ) ' • (3.18) 



Hence 



r^ik ’ 5 x « <3 - 19> 

X(t) being a Wiener process with zero drift and infinitesimal variance 
t. It follows that if 



X(t) £ 8(t) , for {t : 0 s: t £ F(T) [l-F(T) ] _1 ) 



S(t) 

1 - F ( t) 



then 



*; 8(x(t) ) for { t : 0 s: t s: T}. 
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We consider certain special cases results for which are at hand. 
Linear Boundaries . Let B(t) = a + 0 t, Vt ^ 0; a, 3 > 0. Then 
it is known that 

P{X(x) £ a + 3 t, Vt ^ 0} = 1 - e _2a6 (3.20) 

Hence 

S(t) £ a[l-F (t) ] + 6F(t), Vt ^ 0 

_ 2a ft 

with probability 1 - e , and furthermore 

A^(t) £ NF(t) + /N{a[l-F(t)] + 0F(t)}, Vt ^ 0 (3.21) 

—2ct 8 

with approximate probability 1 - e 

Linear Boundaries, Finite Stopping . Let B(t) be as above; then if X 
is a Wiener process and t > 0 is fixed, 

P{X(T f ) < a + gx 1 , V x 1 < t, X(t) £ x < a} = 



$(^-1 _ e 2«g, x-eT-2tt, (3>22) 

/F /F 

where $ is the standard normal distribution; see Cox and Miller [3], 
pp. 220-221. Now use (3.19): 



P{S(t) £ a[l-F(t) ] + 3F(t), V t £ T, S (T) £ y} = 



F (T) 

l-F(T) * 



y(l+T)-0x 


-2ag . 


y(l+T)-3r-2a 


/F J 


“ e 9 


_ /t _ 




(3.23) 



where t = 
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Thus we can represent the probability that the noise component at time 
T is in any Borel set of (-«, a[l-F(T)] + 3F(T)) for any finite 

T. This result can be carried over to an approximation to the distribu- 
tion of A^(t) , subject to boundary avoidance prior to T. 

Similar problems involving two boundaries are also solvable. 

In addition, the distribution of the time required to cross a boundary 
of the form a[l-F(t)] + $F(t) can be derived, although moments are 
not usually expressible in a simple closed form. 

Application of the above material to a problem in sequential 
testing is described in the thesis by Gwinn [5], and an expanded treat- 
ment is in preparation. 



III-2. Diffusion Approximation of (A^(t) ,Q N (t)) . 

In order to establish a diffusion approximation for Q^(t) > 
the number undergoing service, we must treat the bivariate process 
(A^(t) ,Q^(t) ) , focusing on the stochastic elements 

A^t) - Nx(t) 



and 

Q (t) - Ny(t) 

T (t) = — (3.24) 

N S' 

We first proceed as was done earlier. Letting X(t) = f(t)[l-F(t)] ^ 
denote an individual's arrival rate, and assuming exponential service 
times, we write 



PtA^t) = n,Q N (t) = J I Ajg(t) = n,Q N (t) = j} = 1— { A (t) [N-n]+pj }dt+o(dt) 
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p {A^(t) = n+l,Q N (t) = j+ll^Ct) = n,Q N (t) = j} = X (t) [N-n]dt+o(dt) 

P{ A N < t > = n ’ Q N (t) = J -1 I A N ^ t ) = n >Q N (t) = j} = ujdt+o(dt). (3.25) 

other probabilities being negligible. 

Next form the characteristic function expression 

i6 lV t+dt)+i6 2 Q N (t+dt) ie 1 \(t)+ie Q (t) 

E[e ^ 2N ]=E[e ^ 2N {1-X (t) [N-A^t) ] }x 

10,+ie, -i0 

(1-e X Z ) - yQ N (t){l-e Z } + o (dt) (3.26) 

Then substitute from (3.24), divide by dt -> 0, and substitute as in 
(3.8), letting N «>. In order for a limiting ch.f. tj; to exist for 
(S n ,T n ) it turns out that 

x'(t) = X(t) [l-X(t) ] = f(t) 

and (3.27) 

y’(t) = X (t) [ l-x(t) ] - py(t) 

and solution gives the deterministic components of the previous section. 
Finally, ip satisfies 

If = - j[( 0 1 +e 2 ) 2 f(t) + 0 2 yy(t)]* - ( 0 j+6 2 ) X ( t ) g|- Ip - 0 2 P g|- \p (3.28) 

This partial differential equation is the Fourier transformed version 
of a bivariate non-stationary Ornstein-Uhlenbeck process. 

A stochastic differential equation approach runs as follows. 
Approximately, 



dQ N (t) = {X(t) [N-A N (t)] - pQ N (t)}dt + Z(t)/p(t)[N-A N (t)]dt+ Q N (t)dt, 



(3.29) 
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so when (3.24) is applied and N » we find that the deterministic 
part, y, must satisfy 



or 



y ' = f(t) 

t 

» 

y(t) = 



- yy 

e -y(t-s) 



f (s)ds 



(3.30) 



0 

in accordance with earlier results, while the limiting noise is described 
by 



dT = -[X(t)S(t) + pT (t) ] dt + Z(t)/[f (t)+py(t)]dt (3.31) 



Thus T is conditionally Ornstein-Uhlenbeck, given S(t). Also, the 
purely random parts of dT and dS are correlated, as is clear from 
the fact that an additional arrival lengthens the queue; the correlation 
term is f(t), and we can write 

dS (t) = -X (t)S (t)dt + Z (t)/f (t)dt 

dT(t) = -[pT(t)+X(t)S(t)]dt + Z 1 (t)/f(t)dt + Z 2 (t)v^py(t)dt (3.32) 

where and represent independent purely random components, i.e. 

Wiener process derivatives. 



1II-3. Moments. 

Since {S(t),T(t)} is a bivariate diffusion with zero mean we 
know that ip is of the Gaussian form 

^(9 1 ,0 2 ,t) = exp{- j 6' 10) (3.33) 

where 0/ = (0 ,0 2 ) , and is the covariance matrix. Hence, we need 

merely equate coefficients of 0^0^ = 1>2) in order to derive 
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the covariance matrix elements. Not surprisingly, these agree with 
earlier exact results: 



Var [ S (t) ] = x(t) [l-x(t) ] ; Var[T(t)] = y (t) [ 1-y (t) ] 



- F(t) [ 1-F (t) ] 



Cov[S(t) ,T(t) ] = [l-F(t) ] 



-y(t-s) c . N , 
e f (s)ds . 



( 3 . 34 ) 
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IV . Single-Server Models . 

Many of the transitory situations that are encountered in 
practice actually involve competition for facilities; e.g. for repairmen, 
clerks, freeway space, or computers. Therefore in this section we 
describe models for a single-server system confronted by transitory 
demand . 

Once again, a bivariate stochastic process is required to describe 
our model. If N individuals wait to make application for service, 
then the number of new arrivals in (t,t+dt) depends upon the random 
number that has already appeared by time t. The waiting line develop- 
ment is, of course, influenced by the service rate as well. It is the 
coordinated effects of these processes that we represent, and attempt 
to make comprehensible by means of approximations. 

Model 1 . Markovian Arrivals and Exponential Service. 

Let A(t) be the number of arrivals that have occurred by t, 
and Q(t) be the number waiting and in service at t at the service 
facility. Then let 

P. (t) = P{Q(t) = j, A(t) = n | A(0) = 0, Q(0) = 0}, (4.1) 

jn 

for 0 £ j £ n £ N. Given A(t) = n, let X n dt be the probability 
of exactly one arrival in (t,t+dt), and let pdt be the probability 
of exactly one departure in the same interval. The usual independence 
assumptions, cf. Feller [4], are made. Then we can write down the 
forward differential equations for {P^Ct)}: 
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p oo (t) - - x o p oo (t > 



p j n (t) . -<x n+u )p jn (t) + PP J+lin Ct) + x^ 



j = 1,2, .. . ,n-l ; 



P i„ (t) ’ - X „ P 0„ (t) + " P l„<‘> 



P'(t> - -U+|i)P ft) + X P ,(t), 

nn n nn n-1 n-l,n-l 



(4.2) 



where initial conditions, e.g. Pqq(O) = 1, and P (0) = 0 otherwise 
must be specified. The arguments leading to (4.2) are standard and 
will be omitted. Of course, one can also allow arrival and departure 
rates to depend upon t. 

For any explicitly specified X the solution to (4.2) can 

n 

be computed numerically , for (4.2) is only a system of finitely many 
differential equations. If desired, one can Laplace transform through- 
out, but there is no simple explicit formula even for the transform. 

See the paper by Perlas, [11] , for further details. 



Model 1-A . Independent Exponential Arrivals. 

A model of some interest assumes that X = X (N-n) , or, 

n 

equivalently, that each individual’s arrival time is independent with 
distribution 1 - e ^ t . Then a deterministic approximation , that one 
would be tempted to use for N large, is suggested: 

djc(t) = X[N-x(t)]dt (4.3) 

and 

y(t+dt) = max{y_(t) + A[N-x(t)]dt - ydt,0}, (4.4) 
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where x and ^ represent deterministic approximations to A and Q. 
Solutions, subject to A(0) = Q(0) = 0, are 



More will be said about these approximations in the next section. 

Numerical solutions of the above equations have been carried 
out to compare the total waiting times (i) under the full stochastic 
model, (4.1), and (ii) under the deterministic approximation (4.2), (4.3). 
These are summarized in Tables 1 and 2, which show that the deterministic 
approximation improves greatly as N + «. Not surprisingly, the deter- 
ministic approximation falls below the stochastic expectation. 

IV-1. The Probability of No Wait. 



tion focuses upon the probability that no arrival must wait. Other 
applications for this model may well exist, e.g. in military situations. 
Let us suppose that calls occur in accordance with the pure birth process 
of Model 1, and let successive service times {S , n = 1,2,...N} be 



independent but arbitrarily distributed in accordance with F^(x) • If 



2 i(t) = N[l-e ^ t ] 



and 



_y(t) = max{N(l-e ^ t ) - pt,0}. 



(4.4) 



In [6], Kabak describes a transitory problem in which atten- 



n 



T^^ represents the time until the i — arrival occurs, then the 
probability of no wait is seen to be that of the following event 



T (2) >S 1’ T (3) T (2) >S 2’ T (4) T (3) >S 3 




Since times between successive arrivals are independently exponential 
in this Markovian model it is immediate that the probability of no wait is 
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N-l 

p = n 

j=i 



0 



-x t N-l 

e J dF (t) = n F (X ) 
J j-1 J J 



where F is the Laplace-Stieltjes transform of F. 

In the case of Model 1-A this specializes to 



P = 



y + X x y + X 2 ••• y + X N _ X 



© 



N-l 



(¥■ + N-l) S' + N-2) ... (£ + l) 



© rffl 



(N-l)! (N-l) X 



as N -*■ “ 



(4.6) 



(4.7) 



(4.7) 



according to Euler's product definition of the gamma function; see Knopp 
[7]. Of course, if service times are constant rather than exponential 
a simple exact solution is possible: 



N-l 

-- I X(N-j) _ X N (N-l) 



In the event that F is a one-sided stable distribution of order ot, 



0 < a < 1, we have essentially 

. N-l 
r\\ ot 



P = e 



-0“ l 0M)° 

W j-1 



a 1+ot 

-("-I - 



~ e 



1+ot 



(4.9) 



Thus even though the stable laws have right tails long enough to engender 
infinite first moments, the mass near zero governs the above probability. 
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IV-2. The Duration of the Process. 

If we let x denote the time to complete the process of arrival, 
queueing, and service, then it is clear from (4.1) that 

P{t £ t} « P QN (t). (4.10) 

Unfortunately, this distribution has no neat closed form, although 

numerical properties can be derived by solution of (4.2). 

It is worth noting that for large N and under the assumption 

of Model I -A (4.10) is very close to a gamma distribution with mean 
N 

and variance ^r. The explanation is that for large N a queue 
quickly develops behind the server, and does not vanish until the last 
service occurs. The server effectively continues to be active through 
the service of all N arrivals. 
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V . Diffusion Approximation for Single-Server Models . 

The setup to be analyzed now is as follows. N individuals 
independently choose their arrival times from a distribution F(x) , 
the latter possessing a smooth density f(x). In the order of their 
arrival these are served at a single servicing facility; service times 
are exponential at rate Np . Let A^(t) denote the number of arrivals 
that have occurred on or before t, Q^(t) t ^ le number awaiting or 



undergoing service at t, and D^(t) = 



0 



Q^(s)ds is the total delay 



accumulated up to time t. We shall analyze the birth and death process 
(A^(t) ,Q^(t) } ; if F is exponential this is just the process discussed 
in the last section. A useful approximate approach is to derive a 
diffusion approximation, valid for large N. In the present approxima- 
tion, F(t) > pt for a substantial range of t-values (heavy traffic) 
is also a requirement. 



V-l. Deterministic Approximation of (A^Ct) ,Q N (t) > D N (t) } . 

Begin by isolating the deterministic component of the process 
Q^(t) , since that of A^(t) has already been treated. Notice that if 

V° 

— ^ ► y(t) as N -»■ » and it is assumed that y has a derivative, 

then 

Q N (t+h) - Q N (t) = X(t) [N-A N (t)]dt - N U dt + o(dt) (5.1) 

if > 0, while 

Q (t+h) - Q N (t) = X(t)[N-A N (t)]dt + o(dt) (5.2) 

if Q N = 0. Dividing by dt, and by N, and taking the limit as dt -> 0 
and N ->■ 00 yields the differential equation 
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y'(t) = X(t)[l-x(t)] - ph[y(t) ] 



(5.3) 



where y(0) = 0 and 



h(x) = 



1, x > 0 



0, x a 0; 



(5.4) 



the solution is 



y(t) = max{F(t) - pt,0}. 



(5.5) 



It also follows that z(t) , the deterministic component of 
D(t) , is 



z(t) = 



y (s)ds 



(5.6) 



V-2. Diffusion Approximations. 

Next define the supplementary randomness or noise components 
of our three processes: 

A f ,(t) - Nx(t) 

S N (t) & ’ 



Q (t) - Ny(t) 

T (t) = -2 

N ^N 



D (t) - Nz (t) 
UN(t) ^ 



(5.7) 



Begin by writing down the infinitesimal transition probabilities: 



HA^t+dt) = ^(t) , Q N (t+dt) = Q N (t) ,D N (t+dt) = D N (t)+dtQ N (t) | 
A N (t) ,Q N (t) ,D N (t) } = l-X(t) [N-A N (t)]dt-pNh[Q N (t) ]dt+o(dt) , 
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P(Vt+dt) = A N (t)+l,Q N (t+dt) = Q N (t)+l,D N (t+h) = D N (t)+dt[Q N (t)+| | 



A N (t),Q N (t),D N (t)} = X(t)[N-A N (t)]dt+o(dt) 

(5.8) 

?{A^(t+dt) = \( t),Q N (t+dt) = [Q N (t)-l]h[Q N (t)],D N (t+dt) 

= D N (t)+dt[Q N (t)-|] | A N (t),Q N (t),D N (t)} = yNh[Q N (t)]dt+o(dt) 



where once again h(*) represents the unit step at the origin. Now 



put for the joint ch.f. 



® 5 0 2 > ® 3 > E[e 



ie 1 s N (t)+ i0 2 T N (t)+i6 3 U N ( t) 



(5.9) 



In order to derive a partial differential equation for ^ follow the 

conditioning argument used earlier on A^(t) in section III. Also, 

assume throughout that h[Q^(t)] =1, so our approximate solutions will 

not apply when is small (we must avoid the very end of the process, 

and for some distributions, F, the very beginning). Omitting the 

tedious algebra we state that if a limiting noise distribution exists 

its ch.f. ij ) = lim , must satisfy the differential equation 
N-**> N 



■|| = - y[f (t) ( 0 1 +e 2 ) 2 + y 0 2 ] ^ - *(t) (O-j+O^ -jj— 63 g|— (5.10) 



This equation is the Fourier-transformed forward equation of a trivariate 
non-stationary diffusion. As we already know, S(t) is non-s tat ionary 
O.U. Conditionally upon S(t), T(t) is a non-stationary Wiener process. 
This is independently established by deriving the stochastic differential 



equation : 
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dQ N = {X (t) [N-Nx(t)->¥S N (t) - Np}dt 

+ Z(t)/{X(t) [N-Nx(t)-v / NS N (t)+N y dt (5.11) 

= >¥ dT N + Ndy 

After passage to the limit as N + <» and excision of the deterministic 
component the expression 

dT = -X (t)S (t)dt + Z ' (t) /{ f ( t)+p }dt (5.12) 

results. We note that (5.10) may be derived directly from (3.15) and 
(5.12), the stochastic differential equations for the arrival and queue 
noise components. Forgetting about the U component, write 

^(0^,e 2 > t+dt) - ip(e 1 ,e 2 ,t) = 

ie 1 {s(t)+dS}+ie 2 {T(t)+dt) ie 1 s(t)+ie 2 T(t) 

E[e - e ] 

i9 dS+i0 9 dt i0 S(t)+i0 9 T(t) 

= E[{e 1 2 - l}{e 1 }]. 

Now make use of (3.15) and (5.12), and the additional fact that the 
purely random components are correlated : 

-i6 X(t)S(t)dt-i6 2 X(t)S(t)dt 

if)(0 1 ,0 2 ,t+dt)-^(e i ,e 2 ,t) - E[{e 

- \ {f(t)(0 1 +0 9 ) 2 +p0 9 }dt i0 S(t)+i0 9 T(t) 

e 2 12 2 -l}x e 1 2 ] (5.14) 



27 



Next expand exponents to order dt, divide by dt, and let dt -*■ 0; 
after a little rearrangement (5.14) turns into ( 5 . 10 ), with 0 ^ = 0 
and hence the last term removed. Note that the purely random parts of 
S and T are correlated with covariance equal to f(t)dt. Thus one 
might write the equations as follows: 

dS = -A(t)S(t)dt + Z (t)/f (t)dt 

dT = -X (t) S (t) dt + Z x (t)/f (t)dt + Z 2 (t)/yd7 (5.15) 

where Z^ and are independent and purely random. It now follows — 

as it does also by setting 0 ^ = - 6 ^, and 63 = 0 in ( 5 . 10 ) — that 
the noise process T(t) - S(t) is a Wiener process with zero drift and 
infinitesimal variance p. 

In summary, our assumptions (which are akin to those of heavy 
traffic theory in the sense that boundary effects are essentially ignored) 
imply that arrivals, number in system, and total delay are jointly non- 
stationary Gaussian. The actual covariance function will be derived 
subsequently, and some further approximations suggested. From (5.15) 
it is easily seen that a simulation of our approximation can be carried 
out by discretizing time: if t = nA, n = 1,2,..., then 

S ( (n+1) A) = S (nA) - X(nA)S(nA)A + Z^n) /f(nA)A 

T ( (n+1) A) = T(nA) - X(nA)S(nA)A + Z 1 (n)/f(nA)A + Z^n)/^ (5.16) 

where here (Z^n)} and {Z 2 (n)} are mutually independent sequences 
of independent N(0,1) random variables, and S(0) = T(0) = 0. Notice 
too that the derivation of ( 5 . 12 ) may be carried out if the service rate, 
p, is a sufficiently smooth function of t. 
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V-3. Moments. 

The diffusion process characterized by (5.10) is zero mean 
trivariate Gaussian, so 



iKe 1 ,e 2 ,0 3 ,t) = ex P {- \ e'z_Q). 



(5.17) 



where 0/ = [ 0 -^, 02 > 02 -l» — i- ts transpose, and the covariance 

matrix. Now if we differentiate through with respect to t and 0^, 
in accordance with (5.10) and equate coefficients of 0^0^ , we obtain 
the following moments 



Var[S] = F (t) [l-F(t) ] , Var[T] = F(t)[l-F(t)] + yt 



Cov[S,T] = F(t) [l-F(t) ] . (5.18) 



Even more easily, these moments can be obtained from first principles, 
making use of the fact that Var[T-S] = yt as implied by (5.15). The 
fact that the boundary at zero has been neglected, and hence that our 
model is inadequate late in the process, becomes evident in the expres- 
sion for Var[T], which should apparently approach zero with large t, 
but does not. The model should however be accurate near the peak of 
the process. 
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VI . Conclusion . 

The purpose of this paper is to formulate and explore a new class 

of non-s tationary service models. We have shown that under certain 

conditions our models may be described by Gaussian diffusion processes, 

and have hinted at the manner in which known results for such diffusion 

may be utilized to study model properties; for examples we exhibit 

boundaries below which the process moves with prescribed probability, 

and derive equations for the total time spent waiting by all arrivals. 

Although all technical details of approximation accuracy are not yet 

well understood, further theoretical and numerical explorations are 

under way. In any case, the diffusion approximations suggested appear 

to make an offer that the modeler of stochastic phenomena can f t refuse. 

A brief mention of further related questions that deserve attention 

would include (i) the (approximate) distribution of max Q(t), both in 

O^t ^ 00 

sufficient and single server contexts, (ii) approximations when the 
basic process is non-Markovian, (iii) the introduction of decision 
variables, e.g. when to begin servicing, control of arrivals, etc. 

These questions, plus others, are presently under active investigation 



by the authors. 
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Table 1: 


X = 0 


.25 


r— 1 

II 








N 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


E[W] 


7.6 


27.3 


68.7 


136.4 


230.3 


350.0 


494.6 


664.5 


859.4 


1079.4 


W D 


0.1 


13.7 


54. 0 


120.6 


212.7 


330.1 


472.5 


640.0 


832.5 


1050.0 



Table 2: X = 0.75 p = 1 



N 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


E[W] 


10.7 


43.3 


101.5 


184.8 


293.1 


426.4 


584.7 


768.1 


976.4 


1209.7 


W D 


6.0 


36.7 


92.5 


173.3 


279.2 


410.0 


565.8 


746.7 


952.5 


1183.3 



Numerical Comparisons Between E[W] and 



(Deterministic Approximation) 
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